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The Asteroid Redirect Mission (ARM) proposes to retrieve a near-Earth aster¬ 
oid and position it in a lunar distant retrograde orbit (DRO) for later study, crewed 
exploration, and ultimately resource exploitation. During the Caltech Space Chal¬ 
lenge, a recent workshop to design a crewed mission to a captured asteroid in a DRO, 
it became apparent that the asteroid’s low escape velocity (< lcms -1 ) would permit 
the escape of asteroid particles during any meaningful interaction with astronauts or 
robotic probes. This Note finds that up to 5% of escaped asteroid fragments will cross 
Earth-geosynchronous orbits and estimates the risk to satellites from particle escapes 
or complete disruption of a loosely bound rubble pile. 


I. Introduction 

The Asteroid Redirect Mission (ARM) was the result of an initial feasibility study in 2012 by 
the Keck Institute for Space Studies |T], proposing to retrieve a near-Earth asteroid and position 
it in a stable lunar orbit for later study, crewed exploration, and ultimately resource exploitation. 
The ARM is a bold concept, combining both emerging technologies and science, ranging from 
planetary defense to solar system formation to In-Situ Resource Utilization for crewed deep space 
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missions. Requirements derived from expected AP capability and in-orbit stability dictate that a 
captured asteroid be placed in a Lunar distant retrograde orbit (DRO), or possibly Lunar Circulating 
Eccentric Orbits j2J. 

DROs were originally found by Henon [3], who referred to them as the / family, and typically 
encircle the Li and L 2 libration points. Proposed mission concepts exploiting DROs have focused 
on the Earth-Moon |TJ and Jupiter-Europa [5] systems. The family of periodic orbits obtained in 
the circular restricted three body problem (CR3BP) is stable, but the stability region is immersed 
in a region of very unstable motion [B]. Some DROs around the Moon have been found to be stable 
for over 100 years, with the Sun’s gravitational attraction the major external perturbation acting 
on the system [7j . The particular dynamics of the CR3BP has motivated recent studies on transfers 
to lunar DROs using ballistic capture manifold theory mini. 

During the Caltech Space Challenge, a recent workshop to design a crewed mission to a captured 
asteroid in a DRO, it became apparent that the low escape velocity (< lcms -1 ) would permit the 
escape of asteroid particles during any meaningful interaction with astronauts or robotic probes. 
This Note addresses the likely fate of escaped asteroid fragments and estimates the risk to Earth¬ 
orbiting satellites in the event of complete disruption of a loosely bound rubble pile. 

II. Stability of perturbed DROs 

In general, given a certain timescale, DROs possess a region of stability in phase space. Beyond 
this region, a perturbed particle will enter a resonance with the Moon that eventually (months to 
years) destabilizes the orbit. Subsequent motion is chaotic. Particles that escape the DRO enter 
orbits in cis-Lunar space for a period of time until they either collide with the Moon, or are ejected 
from the system entirely, as analysed in Section |III| 

Fig. 0 shows examples of orbits in the CR3BP over a ten year period. A nominal DRO (C = 
2.95) is perturbed by applying in-plane APs when crossing the y-axis close to L\. Four cases are 
distinguished depending on the trajectory of the perturbed particles: a) the DRO remains stable 
over the considered time span, b) the perigee falls below GSOs, c) the particle impacts the Moon 
without crossing GSO, d) the particle escapes the Earth-Moon system without crossing GSO. Note 
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that Lunar flybys are required to enter cis-Lunar regimes. 



x [Earth radii] 


x [Earth radii] 


Fig. 1: Example of orbits in the CR3BP: a) Stable orbit, b) Resonant GSO crossings, c) Moon impact, d) 
Escape without crossing GSO nor impacting the Moon. 


Fig-i shows the stability regions at the DRO’s four cardinal points generated by perturbing 
A Vs on a uniform, symmetric grid. Results of more general DROs were analysed statistically and 
presented in Table[2] Qualitatively, the size of the stability region is similar, but the orientation and 
shape varies with orbital phase. The stability region grows with the Jacobi constant. In addition, 
for smaller values of C the differences in the shape of the stability region across the cardinal points 
are more significant. 
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Fig. 2: Stability maps and fate of particles at the DRO’s cardinal points. A region of stability is 0(10 ms x ) 
across. 


III. Quantifying the risk 

Table [T] studies the effect of the magnitude of the perturbation on the fate of the particles, 
showing the results from a campaign of Monte Carlo simulations. The higher the AV, the more 
particles are ejected from the DRO through a Lunar flyby, increasing the number of GSO crossings 
and Moon impacts. 


4 










Population 

< 10 ms 1 

< 25 ms" 1 

< 44 m s 1 

Population size 

830 

5404 

10000 

Escapes 

10.1% 

15.8% 

22.9% 

GSO 

2.2% 

3.4% 

5.3% 

Moon impact 

10.7% 

19.9% 

30.2% 

Earth impact 

0.0% 

0.0% 

0.0% 

Stable 

77.0% 

60.9% 

41.6% 


Table 1: Percentages of particle fates grouped by maximum perturbing AV on a DRO defined by C = 2.95. 
While smaller AEs result in fewer GSO crossing particles, the cutoff is smaller than 5 ms -1 . 

Table [2] presents the statistical results obtained for simulations based on different DROs, param¬ 
eterized by the Jacobi constant C. Results in the table correspond to the mean value across the four 
cardinal points. We found that the results were robust with respect to applying the perturbation 
at different phases of the orbit, where statistical differences below ±0.5% were observed. This is 
due to the non-chaotic nature of initially perturbed particles. While orbits with a higher Jacobi 
constant are more stable, achieving captured asteroid insertion requires unattainable AV. Persis¬ 
tence time is also important; some of the candidate particles spent months within GSO over a 10 
year period, vastly increasing their interaction cross section compared to transient meteor showers 
or interplanetary dust. 
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c 

(r) [km] 

Period [days] 

GSO-crossing 

Moon impact 

Vi [ms -1 ] 

(N c ) 

N max 

2.86 

156898 

21.59 

6.81% 

5.23% 

-663.93 

41 

864 

2.89 

124372 

18.86 

4.99% 

5.10% 

-590.30 

49 

906 

2.92 

91349 

14.98 

4.02% 

6.40% 

-524.67 

62 

930 

2.95 

65810 

11.10 

5.10% 

13.74% 

-485.44 

66 

859 

2.98 

49458 

8.21 

2.32% 

7.70% 

-473.21 

67 

888 

3.01 

39149 

6.27 

0.96% 

5.18% 

-477.58 

76 

710 

3.04 

32214 

4.94 

0.18% 

1.70% 

-490.29 

83 

876 

3.07 

27290 

4.00 

0.00% 

0.00% 

-507.87 

- 

- 

3.10 

23600 

3.30 

0.00% 

0.00% 

-527.81 

- 

- 


Table 2: Distribution of particle fates for different lunar DR.Os and (AV X , AV y ) £ [—50, 50] x [—50, 50] ms' 1 . 
Each orbit is defined by its mean radius, (r), period, and velocity at the far side of the Moon, i/f. Resonant 
returns are characterized by the average, (N c ) and overall maximum, lV max , of times that particles intersect 
GSO. 


Simulations show that the particles that get to GSO distance from the Earth do not reach low- 
Earth orbit. Risk analyses shall focus on the potential impact of asteroid fragments with satellites in 
GSO. The ARM preliminary mission design considers the early 2020s as the baseline for arriving to 
the lunar DRO. The inclination of the Moon’s orbit with respect to the Earth’s equator determines 
the incoming inclination of the potential impactors, and therefore the satellite orbits most likely to 
be affected. The relative inclination of the Moon reaches a minimum value of 18.2° in Apr-2015, and 
a maximum of 29.0° by Oct-2024. Close approaches detected in the Earth-Moon planar system may 
only interfere with GSO orbits at node crossings. Figure [3] first studies the azimuthal distribution 
of particles at closest approach in the rotating reference. For particles entering GSO there exist 
four differentiated clusters, centered at 0°, ±45° and 180°. We found that the mean time of flight 
to closest approach is roughly 23 months, with a standard deviation of 20 months. The orientation 
of the incoming flux in the fixed frame shows that an isotropic distribution of particles may be 
expected in the first four years. Such a distribution suggests that particles reaching GSO distance 
from Earth may eventually cross GEO at the ascending/descending node. 
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Fig- HI also presents the inclinations of incoming particles at closest approach relative to the 
orbital plane of the Moon, when applying out-of-plane AV Z of up to 100ms _1 . The mean inclination 
is 0.2°, with a standard deviation of 5.1°. This translates into an effective radius of the cross 
section of the particle flux of 3760 km (lex). Hence, 5.7% of GEO will be covered by the incoming 
flux, considering both the ascending and descending nodes and the flux cross section at lex. If 
A14 < 50ms -1 the standard deviation reduces to 3.3°, and the coverage on GEO is 3.6% (lex). 



Min distance to Earth [Earth radii] Time since disruption [yrs] 

Fig. 3: Relative orientation of particles at closest approach for (AV x ,AV y ) G [—50,50] x [—50, 50] ms -1 
(planar case) and (AV y ,AV z ) G [—50,50] x [0,100] ms -1 (out-of-plane motion). 

In the event of complete asteroid disruption into small particles, an ARM design reference 500 T 
asteroid would put 25 T of particles inside GSOs, corresponding to a flux of 34kgs -1 . In contrast to 
resonant DRO escapees, typical interplanetary cosmic dust passes through the GSO region in around 
30 minutes, with a baseline flux of around 100kgs" -m Ipso facto the catastrophic disruption of 
a single 5 m, 500 T asteroid would increase dust risk to satellites by more than 30%. With a mean 
persistence time of around a decade, a flux of only 500kg/day of particles from asteroids in DROs 
will more than double the risk to satellites, necessitating the establishment of good practices for 
asteroid exploration and exploitation in cis-Lunar space. 
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IV. Conclusions 


On average we found that 5% of the particles escaping from a lunar DRO intersect geosyn¬ 
chronous altitude, with no immediate risk to low-Earth orbits. The incoming flux has been found to 
be effectively isotropic and might affect satellites in geostationary orbit (GEO) at node crossings. In 
10 yrs particles have been found to intersect GEO up to 900 times, 63 times on average. A uniform 
flux might be expected during a 40 months period (la). Out-of-plane perturbations yield flux cross 
sections which cover 6% (la) of GEO. These results might serve as a statistical basis for detailed 
analyses of collision risk in particular cases, together with analytical studies on the propagation of 
fragment clouds El- 
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